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We show that a measurement of the reaction energy (\/s) dependence of relative hadron reso- 
nance yields in heavy ion collisions can be used to study the phase structure of the dense strongly 
interacting matter created in these collisions, and investigate the origin of the trends observed in 
the excitation functions of certain soft hadronic observables. We show that presence of chemical 
nonequilibrium in light quark abundance imparts a characteristic signature on the energy depen- 
dence of resonance yields, that differs considerably from what is expected in the equilibrium picture. 
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I. INTRODUCTION 

The exploration of the properties of strongly interact- 
ing, dense quark-gluon matter, and specifically, of the 
equation of state, transport coefficients, degree of equi- 
libration, and phase structure, and the dependence of 
these on the energy and system size is one of the main 
objectives of heavy ion research program. A natural ap- 
proach to these challenges is the study of soft particle 
multiplicities produced in these reactions. This provides 
information about the system properties when these par- 
ticles are created (chemical freeze-out conditions), as well 
as about bulk matter properties (e.g. entropy) which can 
probe deep into the birth history of the fireball. 

Statistical mechanics techniques have in this context a 
long and illustrious history [H, 0, SQ. The systematic 
and quantitative comparison of data to the statistical 
hadronization (SH) model is, however, a comparatively 
recent field [!, i, 0, 1, i, O, EH- A consensus has de- 
veloped that the SH model can indeed fit most, if not all 
particle yields measured at experiments conducted at a 
wide range of energies. Measurements conducted at the 
GSI Schwerionen Synchrotron (SIS), BNL's Alternating 
Gradient Synchrotron (AGS),CERN's Super Proton Syn- 
chrotron (SPS),and BNL's Relativistic Heavy Ion Col- 
lider (RHIC) have successfully been analyzed using SH 
ansatze. 

When this consensus is considered more carefully, we 
see that in technical detail the applied SH models dif- 
fer regarding the chemical equilibration condition that 
is presumed. As a result, it has not as yet been possi- 
ble to agree statistical physics, if any, is responsible for 
the striking trends observed in the energy dependence of 
some observed hadronic yields. 

In this paper we will indicate that further progress can 
be made with help of hadron resonances. Hadron reso- 
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nances, such as e.g. A ++ excitation of the proton p differ 
typically from the "stable" particle by internal structure, 
rather than chemical quark content. Hence within the 
SH approach their production is mostly controlled by the 
"temperature" T at which they are created. 



II. EQUILIBRIUM AND NON-EQUILIBRIUM 
FREEZE-OUT CONDITION 



In the SH model there are two types of chemical equi- 
librium [T^|.: all models assume relative chemical equi- 
librium, but some also assume absolute chemical equilib- 
rium which implies the presence of just the right abun- 
dances of valance up, down, and even strange quark pairs. 
There are qualitative differences in the results obtained 
in the description of hadron production with or without 
using the hypothesis of absolute chemical equilibrium: if 
the system of produced hadrons is considered to be in 
absolute chemical equilibrium, then at highest heavy ion 
reaction energy one obtains chemical frcczc-out temper- 
ature T ~ 160 — 170 MeV. Values as low as T ~ 50 MeV 
are reported at lowest reaction energies available. 

The energy dependence of the freeze-out temperature 
than follows the trend indicated in panel (a) of figure 
[1] as the collision energy increases, the freeze-out tem- 
perature increases and the baryonic density (here bary- 
onic chemical potential hb) decreases @. An increase of 
freeze-out temperature with ^fs is expected on general 
grounds, since with increasing reaction energy a greater 
fraction of the energy is carried by mesons created in the 
collision, rather than pre-existing baryons [l9j |. 

Further refinements in the approach described above 
are often implemented and could be of relevance: 

• Allowance for strangeness chemical nonequilibrium 
is necessary to obtain a good description of strange 
particle yields [H, [H, 111 at low V^- Tnis is 
accomplished by introducing strangeness phase 
space occupancy 7 S ; 
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FIG. 1: (Color online) Dependence of freeze-out temperature T and baryo-chemical potential /xb on reaction energy in the 
Equilibrium (panel (a), ||) and non-equilibrium (panel (b),[j|) freeze-out models. The direction of the arrow corresponds to 
increasing ^fs. The equilibrium dependence of T and \xb in the panel (a) is not significantly altered by the introduction of the 
fitted phase space occupancy 7 S and/or the implementation of the Canonical ensemble for strangeness. The "star" in panel (b) 
corresponds to the point where the transition to the supercooled regime occurs and the phase space changes from chemically 
under-saturated (jj q < 1) to chemically over-saturated (7, > 1). This point also corresponds to the energy of the "kink" and 
the tip of the "horn" 



• At smaller reaction energies and for smaller reac- 
tion system sizes, it is likely that the fireball is 
well away from the thermodynamic limit. In this 
case canonical treatment of strangeness is applied 




These effects do not materially alter the behaviour of 
temperature and chemical potential shown in the panel 
(a) of Fig. [Q 

What is most striking in these results is that there is 
no sign of any structure when the reaction energy varies. 
However, there are non-continuous features in the energy 
dependence of hadronic observables, such as the "kink" 
in the multiplicit y p er number of participants and the 
"horn" [{J [20, I2H |22|| in certain particle yield ratios. An 
effort was made to interpret this in terms of a shift from 
baryon to meson dominance [22| of the hadron yields. 
However no matter how hard one tries, in the chem- 
ical equilibrium model even the simple observable like 
K + /ir + remains a smooth function of reaction energy, in 
contrast to the experimental results. Introduction of 7,,. 
and deviations from the thermodynamic limit, while they 
help in bringing some of the model predictions closer to 
the data, has so far not managed to reproduce the sharp- 
ness of features such as the kink and the horn. 

Non-monotonic behaviour of particle yield ratios could 
indicate a novel reaction mechanism, e.g. onset of the de- 
confinement phase [21] . In such a situation, the smooth- 
ness of the chemical freeze-out temperature dependence 
on energy would be surprising, since it would imply that 



at all energies, from about 1 A GeV at SIS, to the high- 
est RHIC values, there is no change in either the fireball 
evolution dynamics, nor any other imprint from the de- 
confined phase on the freeze-out condition, which, how- 
ever is visible in the strangeness and entropy yield that 
K + and, respectively, 7r + represent. 

Furthermore, we note that the fireball of hadronic mat- 
ter formed is a relatively small system, expanding rapidly, 
with its content undergoing a phase transformation, or 
even phase transition. In this complex and rapidly evolv- 
ing circumstance, one could imagine that the absolute 
chemical equilibrium, not always, or even ever, holds. In 
particular, if the expanding system undergoes a fast con- 
version from a Quark Gluon Plasma (QGP) to hadrons, 
chemical non-equilibrium |5;] , and super-cooling [23L |24| 
go hand in hand, due to entropy and flavor conservation 
requirements. 

One can look at this situation again removing the hy- 
pothesis of absolute chemical equilibrium among hadrons 
produced. The systematic behaviour of T with energy in 
this case is quite different [9( , as is shown in panel (b) of 
figure[T] The two higher T values at right are for 20 (low- 
est SPS) and (most to right) 11.6 A GeV (highest AGS) 
reactions. In these two cases the source of particles is a 
hot chemically under-saturated (T ~ 170 MeV ) fireball. 
Such a system could be a conventional hadron gas fireball 
that had not the time to chemicall y e quilibrate. Other 
options were considered in Refs. d, HH, such as a phase 
of constituent massive quarks. 
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Following the thick arrow in panel (b) of figure [T] we 
note that somewhat smaller temperatures are found with 
further increasing heavy ion reaction energies. Here it is 
possible [H, HUto match the entropy of the emerging 
hadrons with that of a system of nearly massless partons 
when one considers supercooling to T ~ 140 MeV, while 
both light and strange quark phase space in the hadron 
stage acquire significant over-saturation with the phase 
space occupancy 7 g= „ !( 2 > 1 and at higher energy also 
7 S > 1. A drastic change in the non-equilibrium con- 
dition occurs near 30 A GeV, corresponding to the dip 
point on right in panel (b) of the figure Q] (marked by 
an asterisk). At heavy ion reaction energy below (i.e. to 
right in panel (b) of figure [1]) of this point, hadrons have 
not reached chemical equilibrium, while at this point, as 
well as, at heavy ion reaction energy above (i.e. at and to 
left in panel (b) of figure [T]), hadrons emerge from a much 
denser and chemically more saturated system, as would 
be expected were QGP formed at and above 30 A GeV. 
This is also the heavy ion reaction energy corresponding 
to the "kink", which tracks the QGP's entropy density 
(higher w.r.t. a hadron gas), and the peak of the "horn" 
[201 ] , which tracks the strangeness over entropy ratio (also 
higher w.r.t. a hadron gas). 

Concluding this discussion, comparing panel (a) with 
(b) we see in figure [1] a quite different behaviour. On 
panel (a), for the chemical equilibrium model, with in- 
crease of the collision energy following the black arrow, 
we see monotonic increase of the chemical freeze-out tem- 
perature, with no hint of new physics in a wide range of 
heavy ion collision energies spanning the range of SIS, 
AGS, SPS and RHIC. On panel (b), we see that when 
relative and absolute chemical equilibrium is considered 
[IH , with yields of individual hadrons satisfying the rel- 
ative, but not the absolute chemical equilibrium, the ex- 
perimental particle yield data is best described with a 
temperature profile as function of reaction energy which 
is not monotonic. There is a minimum value of T, at 
the point when the rapid change of the chemical compo- 
sition of produced hadrons is occurring. This is clearly 
suggestive of a change in the reaction mechanism. 

The main reason for the wider acceptance of the equi- 
librium approach 7, = 1 is its greater simplicity, there 
are fewer parameters. Moreover, considering the qual- 
ity of the data the non-equilibrium parameter 7^ is not 
necessary to pull the statistical significance above it's 
generally accepted minimal value of 5 %. On the other 
hand, the parameters 7 ? and 7 S were introduced on phys- 
ical grounds [HI, [T^, [23| , thus these are not arbitrary fit 
parameters. Moreover, these parameters, when used in 
a statistical hadronization fit, converge to theoretically 
motivated values. They also help to explain the trends 
observed in the energy dependence of hadronic observ- 
ables. Finally, 7 9 > 1 in AA reactions describes the en- 
hancement of the baryon to meson ratio yield at RHIC, 
as compared to elementary interactions, which dynami- 
cally arises in the recombination hadronization at fixed 
hadronization temperature. 



III. RESONANCE RATIOS AS CHEMICAL 
FREEZE-OUT TEMPERATURE PROBES 

In this paper, we propose the energy dependence of 
the resonance yields as a possible experimental observ- 
able, capable to discriminate the two scenarios, chemical 
equilibrium and non-equilibrium, and thus to establish 
the need to use 7^ in statistical hadronization analysis of 
experimental data. 

Many strong interaction resonances, a set we denote 
by the collective symbol R* (such as IT" (892), A(1232), 
E*(1385), A*(1520), 5* (1530) (2f|) carry the same 
valance quark content as their ground-state counter-parts 
R (corresponding: K, N, S, A, 3). R* typically decay 
by emission of a pion, R* — > R + tt. Considering the par- 
ticle yield ratio R*/R in the Boltzmann approximation 
(appropriate for the particles considered), we see that 
all chemical conditions and parameters (equilibrium and 
non-equilibrium) cancel out, and the ratio of yields be- 
tween the resonance and it's ground state is a function of 
the masses, and the freeze-out temperature, with second 
order effects coming from the cascading decays of other, 
more massive resonances [H, [HJ : 

N R . _ 9R-W(2f-) + b jR . g 3 W (g) 

Nr~ g R W +Ek^ R hR9kW U 

where W[x) = x^K^ix) is the (relativistic) reduced one 
particle phase space, K^{x) being a Bessel function, g is 
the quantum degeneracy, and bjR is the branching ratio 
of resonance j decaying into R. 

When we study the results arising from Eq. (fTJ) , we con- 
sider only strong decay contributions, weak decay feed- 
down, such as A — * p, £ — * p, 3 — * A, and f2 — > 3 
has to be eliminated from the data sample. Given that 
existing SPS [27| and RHIC [28| experiments, as well 
as the planned Compressed Baryonic Matter (CBM) ex- 
periment [29l | have both a tracking resolution permitting 
precise primary vertex cuts (weak decay tracks originate 
from points well away from the primary vertex), as well 
as a momentum resolution capable of identifying reso- 
nances [3(J HH 32], this requirement is realistic. 

Because of the radically different energy dependence 
of freeze-out temperature in the scenarios of @ and @, 
seen in figure [TJ the prediction for the resonance ratios 
Eq. vary greatly between these two scenarios. In the 
equilibrium scenario the temperature goes up with heavy 
ion reaction energy, and thus the resonance abundance 
should go smoothly up for all resonances. On the other 
hand, the nonequilibrium scenario, with a low tempera- 
ture arising only in some limited reaction energy domain, 
will lead to resonance abundance which should have a 
clear dip at that point, but otherwise remain relatively 
large. 

We have evaluated several resonance relative ratios 
shown in figure [5] within the two scenarios, using the sta- 
tistical hadronization code SHARE [H, [34| • For the non- 
equilibrium scenario, we have used the parameters given 
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in [9(, table I. For the equilibrium scenario, we used the 
parametrization given in Q figures 3 and 4. In the latter 
case, the strangeness and isospin chemical potentials [i s 
and [ii3 were obtained by requiring that net strangeness 
be zero, and net charge per baryon in all particles pro- 
duced be the same as in the colliding system. We have 
performed spot checks of the validity of the statistical 
parameters used and found that under the assumptions 
made these are the best parameter sets. 

As seen in figure [U the expected trend with a/s is ap- 
parent in all considered resonance ratios, though in cases 
where the mass difference is large, the effect is much more 
pronounced than in some others. Indeed, for many of the 
ratios we present the experimental error may limit the 
usefulness of our results, however, the opposite energy 
dependence may prove to be helpful in discriminating 
the behaviour. In principle we have presented 12 differ- 
ent ratios, though in some cases the difference between 
particle and antiparticle ratio is small. It is smallest in 
cases when the baryochemical potential has neither a di- 
rect, nor a large indirect influence, and the difference in 
light quark flavor content is the smallest. 

It is interesting to note here that the often ignored 
quark flavor effect (isospin effect) is responsible for most 
of the difference between particle and antiparticle ratios. 
This is at first a counter intuitive result, but it can be 
understood in quantitative simple manner. We recall the 
relation u/d oc A/3 cx d/u where u,d refer to the yield 
of up, and, respectively, down valance quarks. In cases, 
such as e.g. E* — > A the u, d valance quark content is 
different for the R* and R particles, leading to A/3 de- 
pendence of the R* /R ratio. High mass cascading res- 
onances, where the strangeness content can be different 
(e.g. S(1690) — > KH), are a further source of difference 
between particle and antiparticle resonance ratios, espe- 
cially so in the high baryo-chemical potential regime. 

Many of the experimental data points needed in the 
analysis presented in 0] and @ are to this day still pre- 
liminary. This means that some of the results we rely 
on could in the end be somewhat different. However, be- 
cause of the cancellation (to a good approximation) of 
the baryo- and strangeness chemical potentials, the qual- 
itative prediction for the energy dependence of the reso- 
nance yields within the two models is robust. Namely, 
within the chemical equilibrium model the temperature 
of chemical freeze-out must steadily increase and so does 
the R* I R ratio. For the chemical non-equilibrium model 
the R* I R dip primarily relies on the response of T to 
the degree of chemical equilibration: prior to chemical 
equilibrium for the valance quark abundance, at a rela- 
tively low reaction energy, the freeze-out temperature T 
is relatively high. At a critical energy, T drops as the 
hadron yields move to or even exceed light quark chemi- 
cal equilibrium, yet reaction energy is still not too large, 
and thus the baryon density is high and meson yield low. 
As reaction energy increases further, T increases and the 
R* /R yield from that point on increases. The drop in 
R* j R at critical T, would be completely counter-intuitive 



in an equilibrium picture. It would hence provide over- 
whelming evidence that non-equilibrium effects such as 
supercooling, where such a drop would be expected, are 
at play. We further argue that such a drop can not be 
reproduced by resonance rescattering/regeneration. 

Pseudo-elastic processes such as Rir — > R* — > Rtt and 
post-decay R* — ► Rn scattering of decay products in mat- 
ter could potentially considerably alter the observable fi- 
nal ratio of detectable R* to R. The combined effect of 
rescattering and regeneration has not been well under- 
stood. We have argued that the formation of additional 
detectable resonances is negligible [35j | , while scattering of 
decay products can decrease the visible resonance yields 
except for sudden hadronization case. Other groups have 
studied this in quantitative manner. 

Assuming a long lived hadron phase, the energy de- 
pendence of most of the resonance ratios considered here 
has been calculated in a hadronic quantum molecular 
dynamics model. The result (figure 7 and 8 in [36|]) is 
qualitatively similar to the chemical equilibrium results 
for resonance ratios, we see a smooth rise with energy. 
Thus, in the case of chemical equilibrium, with a consid- 
erable separation between chemical and thermal freeze- 
out inherent in Ref. [3t| rescattering and regeneration 
will affect the quantitative R* / R ratio, but will not al- 
ter the dependence on heavy ion reaction energy shown 
in figure [5] On the other hand, chemical nonequilibrium 
implies absence of a long lived hadron phase. Because of 
this, the calculation [36] would not be applicable and the 
resonance abundance should be in closer quantitative, as 
well as qualitative agreement, with the predictions of fig- 
ure [2J Rescattering and regeneration, therefore, should 
not alter the predicted pattern of either the equilibrium 
model, where all R*/R should rise with energy, or the 
non-equilibrium model, where at the critical energy all 
R* j R should experience a dip. 

If both of these predictions prove inaccurate, and the 
R* j R abundance turns out to be resonance specific with 
no uniform rises and dips as function of energy, this would 
signify that freeze-out is determined by reaction kinetics 
rather than thermodynamic conditions... In this case, R* 
abundance is determined more by Tr*t, where r is the 
lifetime of the fireball, than by ma- /T. 



IV. CONCLUSIONS 

We have argued that a measurement of the en- 
ergy dependence of ratios such as K*/K, A/p, S*/A, 
A(1520)/A, S*/S and other such ratios might be instru- 
mental in clarifying the freeze-out conditions in heavy ion 
collisions, especially at low reaction energy. A resonance 
abundance monotonically rising with energy from the 
AGS energy range would suggest that the best statistical 
description of heavy ion data is based on chemical equi- 
librium, and that as collision energy increases, freeze-out 
temperature rises monotonically. If, on the other hand, 
resonance abundance shows a consistent dip, possibly at 
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FIG. 2: (Color online) Ratio of resonance to the stable particle. Thick lines for particles with strange quark content, thin lines 
for particles with anti-strange quark content, as a function of energy. Solid black lines refer to the equilibrium fits (-y q ,s = 1), 
with the parameters for AGS and SPS energies taken from jjj. Dashed red lines refer to non-equilibrium fits (7q, s fitted), with 
the best fit parameters for AGS and SPS energies taken from [tj . 
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the energy coinciding with the other non-monotonic fea- 
tures recently observed in particle yields and ratios [20j , 
it would be a strong evidence that what we are seeing 
is, at and above this dip, a freeze-out from a supercooled 
high entropy density phase. GT would like to thank C. 
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